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Abstract 

Primordial black holes (PBHs) are a profound signature of primordial 
cosmological structures and provide a theoretical tool to study nontriv- 
ial physics of the early Universe. The mechanisms of PBH formation 
are discussed and observational constraints on the PBH spectrum, or ef- 
fects of PBH evaporation, are shown to restrict a wide range of particle 
physics models, predicting an enhancement of the ultraviolet part of the 
spectrum of density perturbations, early dust-like stages, first order phase 
transitions and stages of superheavy metastable particle dominance in the 
early Universe. The mechanism of closed wall contraction can lead, in the 
inflationary Universe, to a new approach to galaxy formation, involving 
primordial clouds of massive BHs created around the intermediate mass 
or supermassive BH and playing the role of galactic seeds. 
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1 Introduction 

The convergence of the frontiers of our knowledge in micro- and macro- worlds 
leads to the wrong circle of problems, illustrated by the mystical Uhroboros (self- 
eating-snake) . The Uhroboros puzzle may be formulated as follows: The theory 
of the Universe is based on the predictions of particle theory, that need cosmology 
for their test. Cosnioparticle physics [TJ [U [31 HI [SJ [S] offers the way out of this 
wrong circle. It studies the fundamental basis and mutual relationship between 
micro-and macro-worlds in the proper combination of physical, astrophysical 
and cosmological signatures. Some aspects of this relationship, which arise in 
the astrophysical problem of Primordial Black Holes (PBH) is the subject of 
this review. 

In particle theory Noether's theorem relates the exact symmetry to con- 
servation of respective charge. Extensions of the standard model imply new 
symmetries and new particle states. The respective symmetry breaking induces 
new fundamental physical scales in particle theory. If the symmetry is strict, 
its existence implies new conserved charge. The lightest particle, bearing this 
charge, is stable. It gives rise to the deep relationship between dark matter 
candidates and particle symmetry beyond the Standard model. 

The mechanism of spontaneous breaking of particle symmetry also has cos- 
mological impact. Heating of condensed matter leads to restoration of its sym- 
metry. When the heated matter cools down, phase transition to the phase of bro- 
ken symmetry takes place. In the course of the phase transitions, corresponding 
to given type of symmetry breaking, topological defects can form. One can di- 
rectly observe formation of such defects in liquid crystals or in supcrfluid He. In 
the same manner the mechanism of spontaneous breaking of particle symmetry 
implies restoration of the underlying symmetry. When temperature decreases 
in the course of cosmological expansion, transitions to the phase of broken sym- 
metry can lead, depending on the symmetry breaking pattern, to formation of 
topological defects in very early Universe. Defects can represent new forms of 
stable particles (as it is in the case of magnetic monopoles [71 [51 [TUl [TTl [T^ 1 . 
or extended structures, such as cosmic strings [131 [H] or cosmic walls [TS]. 

In the old Big bang scenario cosmological expansion and its initial conditions 
were given a priori [THl HZ] . In the modern cosmology expansion of Universe 
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and its initial conditions are related to inflation [TSl [IHl UHl HIl 112] ) baryosyn- 
thesis and nonbaryonic dark matter (see review in [23l[24]). Physics, underlying 
inflation, baryosynthesis and dark matter, is referred to extensions of the stan- 
dard model, and variety of such extensions makes the whole picture in general 
ambiguous. However, in a framework of each particular physical realization 
of inflationary model with baryosynthesis and dark matter the corresponding 
model dependent cosmological scenario can be specified in all details. In such 
scenario main stages of cosmological evolution, structure and physical content 
of the Universe reflect structure of the underlying physical model. The lat- 
ter should include with necessity the standard model, describing properties of 
baryonic matter, and its extensions, responsible for inflation, baryosynthesis 
and dark matter. In no case cosmological impact of such extensions is reduced 
to reproduction of these three phenomena only. A nontrivial path of cosmolog- 
ical evolution, specific for each particular realization of inflational model with 
baryosynthesis and nonbaryonic dark matter, always contains some additional 
model dependent cosmologically viable predictions, which can be confronted 
with astrophysical data. Here we concentrate on Primordial Black Holes as 
profound signature of such phenomena. 

It was probably Pierre-Simon Laplace [25j in the beginning of XIX century, 
who noted first that in very massive stars escape velocity can exceed the speed 
of light and light can not come from such stars. This conclusion made in the 
framework of Newton mechanics and Newton corpuscular theory of light has 
further transformed into the notion of " black hole" in the framework of general 
relativity and electromagnetic theory. Any object of mass M can become a 
black hole, being put within its gravitational radius rg = IGMjt? . At present 
time black holes (BH) can be created only by a gravitational collapse of com- 
pact objects with mass more than about three Solar mass [SHIHT]. It can be a 
natural end of massive stars or can result from evolution of dense stellar clus- 
ters. However in the early Universe there were no limits on the mass of BH. 
Ya.B. Zeldovich and I.D. Novikov [5S] noticed that if cosmological expansion 
stops in some region, black hole can be formed in this region within the cosmo- 
logical horizon. It corresponds to strong deviation from general expansion and 
reflects strong inhomogeneity in the early Universe. There are several mecha- 
nisms for such strong inhomogeneity and we'll trace their links to cosmological 
consequences of particle theory. 

Primordial Black Holes (PBHs) are a very sensitive cosmological probe for 
physics phenomena occurring in the early Universe. They could be formed by 
many different mechanisms, e.^., initial density inhomogeneities |291 130] and 
non-linear metric perturbations [311 [32l [33] , blue spectra of density fluctuations 
[31 [33 [31 [371 [3H1 [3S|, a softening of the equation of state [H [Ml [35], de- 
velopment of gravitational instability on early dust-like stages of dominance of 
supermassive particles and scalar fields [HI [HI [131 [H] and evolution of gravita- 
tionally bound objects formed at these stages [3S1[3S], collapse of cosmic strings 
da [H [IS [501151] and necklaces [52], a double infiation scenario [53l [M [SSI [56] , 
first order phase transitions [SZl [SS [SS [60] [61] , a step in the power spectrum 
[621 [eg , etc. (see [35l E [1 [Ml [65] for a review) . 
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Being formed, PBHs should retain in the Universe and, if survive to the 
present time, represent a specific form of dark matter [66l [67l [3l l4l [68l [69l [70l [65l 
[71] . Effect of PBH evaporation by S.W. Hawking [72] makes evaporating PBHs 
a source of fluxes of products of evaporation, particularly of 7 radiation [73] . 
MiniPBHs with mass below 10^^ g evaporate completely and do not survive to 
the present time. However, effect of their evaporation should cause influence on 
physical processes in the early Universe, thus providing a test for their existence 
by methods of cosmoarcheology [71], studying cosmological imprints of new 
physics in astrophysical data. In a wide range of parameters the predicted 
effect of PBHs contradicts the data and it puts restrictions on mechanism of 
PBH formation and the underlying physics of very early Universe. On the other 
hand, at some fixed values of parameters, PBHs or effects of their evaporation 
can provide a nontrivial solution for astrophysical problems. 

Various aspects of PBH physics, mechanisms of their formation, evolution 
and effects are discussed in|75l[76l[77l|78l|35l[7a[80l[8ll[8a[83l[8a 
IB8, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99^, ITUni [TOTl [TM Um UMl [TU51 [TM 
QOI 008.109.110. Ill, 112, 113. .114„ 115.1lBl[TT71[TT51[TTOl[T^[m[m[n5] 
particularly specifying PBH formation and effects in braneworld cosmology 
[m [Hil [HI [113, on inflationary preheating [HH], formation of PBHs in 
QCD phase transition |129lll30j . properties of superhorizon BHs ' 13111132) . role 
of PBHs in baryosynthesis [HI [HI [I3S1 [Hi [133 , effects of PBH evaporation 
in the early Universe and in modern cosmic ray, neutrino and gamma fluxes 
[1381 fT39l fT40l nm [T421 [T43l fT44l [T45l [T46l iMTl [T48l fT49l [TSOl [TsT, T^, HSl 
[I541[T551[I561II571[I581[I591[I601[I6I1 168, 169j, in 

creation of hypothetical particles [1701 11711 [172, 173] , PBH clustering and cre- 
ation of supermassive BHs [1741 11751 11761 [T77> .178. .179j , effects in cosmic rays 
and colliders from PBHs in low scale gravity models |1801 1181] . Here we outline 
the role of PBHs as a link in cosmoarcheoLOGICAL chain, connecting cosmo- 
logical predictions of particle theory with observational data. We discuss the 
way, in which spectrum of PBHs reflects properties of superheavy metastable 
particles and of phase transitions on inflationary and post-inflationary stages. 
We briefly review possible cosmological reflections of particle physics (section[5]), 
illustrate in section[3]some mechanisms of PBH formation on stage of dominance 
of superheavy particles and fields fsubsection l3.ip and from second order phase 
transition on inflationary stage. Effective mechanism of BH formation during 
bubble nucleation provides a sensitive tool to probe existence of cosmological 
first order phase transitions by PBHs (section [4]) . Existence of stable remnants 
of PBH evaporation can strongly increase the sensitivity of such probe and we 
demonstrate this possibility in section [5] on an example of gravitino production 
in PBH evaporation. Being formed within cosmological horizon, PBHs seem 
to have masses much less than the mass of stars, constrained by small size of 
horizon in very early Universe. However, if phase transition takes place on infla- 
tionary stage, closed walls of practically any size can be formed fsubsection l6.2p 
and their successive collapse can give rise to clouds of massive black holes, which 
can play the role of seeds for galaxies (section [6]). The impact of constraints and 
cosmological scenarios, involving primordial black holes, is briefly discussed in 
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section [T] 



2 PBHs as cosmological reflection of new physics 

The simplest primordial form of new physics is a gas of new stable massive 
particles, originated from early Universe. For particles with mass m, at high 
temperature T > m the equilibrium condition, n-av-t > 1 is valid, if their anni- 
hilation cross section a > l/{mm.pi) is sufficiently large to establish equilibrium. 
At T < m such particles go out of equilibrium and their relative concentration 
freezes out. Weakly interacting species decouple from plasma and radiation at 
T > m, when n-av-t 1, i.e. at Tdec ~ (cmpj)"^. This is the main idea 
of calculation of primordial abundance for WIMP-like dark matter candidates 
(see e.g. [3J|31[73] for details). The maximal temperature, which is reached in 
inflationary Universe, is the reheating temperature, T^, after inflation. So, very 
weakly interacting particles with annihilation cross section a < l/{Tj.mpi), as 
well as very heavy particles with mass m ^ Tr can not be in thermal equilib- 
rium, and the detailed mechanism of their production should be considered to 
calculate their primordial abundance. 

Decaying particles with lifetime r, exceeding the age of the Universe, tu, t > 
tu, can be treated as stable. By definition, primordial stable particles survive 
to the present time and should be present in the modern Universe. The net 
effect of their existence is given by their contribution into the total cosmological 
density. They can dominate in the total density being the dominant form of 
cosmological dark matter, or they can represent its subdominant fraction. In 
the latter case more detailed analysis of their distribution in space, of their 
condensation in galaxies, of their capture by stars. Sun and Earth, as well 
as effects of their interaction with matter and of their annihilation provides 
more sensitive probes for their existence. In particular, hypothetical stable 
neutrinos of 4th generation with mass about 50 GeV are predicted to form 
the subdominant form of modern dark matter, contributing less than 0,1 % to 
the total density [1861 1187j . However, direct experimental search for cosmic 
fluxes of weakly interacting massive particles (WIMPs) may be sensitive to 
existence of such component (see |182[ 11831 1184[ 1185] and references therein). 
It was shown in [188[ 11891 11901 I191j that annihilation of 4th neutrinos and 
their antineutrinos in the Galaxy can explain the galactic gamma-background, 
measured by EGRET in the range above 1 GeV, and that it can give some clue 
to explanation of cosmic positron anomaly, claimed to be found by HEAT. 4th 
neutrino annihilation inside the Earth should lead to the flux of underground 
monochromatic neutrinos of known types, which can be traced in the analysis 
of the already existing and future data of underground neutrino detectors [190|. 
[T92l[T93l[T94] . 

New particles with electric charge and/or strong interaction can form anoma- 
lous atoms and contain in the ordinary matter as anomalous isotopes. For 
example, if the lightest quark of 4th generation is stable, it can form stable 
charged hadrons, serving as nuclei of anomalous atoms of e.g. anomalous he- 
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Primordial unstable particles with lifetime, less than the age of the Universe, 
T < tu, can not survive to the present time. But, if their lifetime is sufficiently 
large to satisfy the condition r 3> (mpi/m) ■ (1/m), their existence in early Uni- 
verse can lead to direct or indirect traces. Cosmological flux of decay products 
contributing into the cosmic and gamma ray backgrounds represents the direct 
trace of unstable particles. If the decay products do not survive to the present 
time their interaction with matter and radiation can cause indirect trace in the 
light element abundance or in the fluctuations of thermal radiation. 

If particle lifetime is much less than 1 s multi-step indirect traces are possible, 
provided that particles dominate in the Universe before their decay. On dust-like 
stage of their dominance black hole formation takes place, and spectrum of such 
primordial black holes traces particle properties (mass, frozen concentration, 
lifetime) [35] . Particle decay in the end of dust like stage influences the baryon 
asymmetry of the Universe. In any way cosmophenomenoLOGICAL chains link 
the predicted properties of even unstable new particles to the effects accessible 
in astronomical observations. Such effects may be important in analysis of the 
observational data. 

Parameters of new stable and metastable particles are also determined by a 
pattern of particle symmetry breaking. This pattern is reflected in a succession 
of phase transitions in the early Universe. First order phase transitions proceed 
through bubble nucleation, which can result in black hole formation (see e.g. 
[59] and [65] for review and references). Phase transitions of the second order 
can lead to formation of topological defects, such as walls, string or monopoles. 
The observational data put severe constraints on magnetic monopole [9] and cos- 
mic wall production [T^ , as well as on the parameters of cosmic strings [13l [14] . 
Structure of cosmological defects can be changed in succession of phase tran- 
sitions. More complicated forms like walls-surrounded-by-strings can appear. 
Such structures can be unstable, but their existence can leave a trace in non- 
homogeneous distribution of dark matter and give rise to large scale structures 
of nonhomogeneous dark matter like archioles [2011 12021 1203j . Primordial Black 
Holes represent a profound signature of such structures. 

3 PBHs from early dust-like stages 

A possibility to form a black hole is highly improbable in homogeneous ex- 
panding Universe, since it implies metric fluctuations of order 1. For metric 
fluctuations distributed according to Gaussian law with dispersion 



a probability for fluctuation of order 1 is determined by exponentially small tail 
of high amplitude part of this distribution. This probability can be even more 
suppressed in a case of non-Gaussian flutuations [21] . 
In the Universe with equation of state 




(1) 



P = le, 



(2) 
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with numerical factor 7 being in the range 

< 7 < 1 (3) 

a probabihty to form black hole from fluctuation within cosmological horizon is 
given by (see e.g. [31 U] for review and references) 

WpBHCcexp(^-^L_y (4) 

It provides exponential sensitivity of PBH spectrum to softening of equation of 
state in early Universe (7 — )■ 0) or to increase of ultraviolet part of spectrum of 
density fluctuations (((5^) — > 1). These phenomena can appear as cosmological 
consequence of particle theory. 



3.1 Dominance of superheavy particles in early Universe 

Superheavy particles can not be studied at accelerators directly. If they are sta- 
ble, their existence can be probed by cosmological tests, but there is no direct 
link between astrophysical data and existence of superheavy metastable parti- 
cles with lifetime r <C Is. It was first noticed in [41] that dominance of such 
particles in the Universe before their decay at t < t can result in formation of 
PBHs, retaining in Universe after the particles decay and keeping some infor- 
mation on particle properties in their spectrum. It provided though indirect 
but still a possibility to probe existence of such particles in astrophysical obser- 
vations. Even the absence of observational evidences for PBHs is important. It 
puts restrictions on allowed properties of superheavy metastable particles, which 
might form such PBHs on a stage of particle dominance, and thus constrains 
parameters of models, predicting these particles. 
After reheating, at 

T <Ta=rm (5) 

particles with mass m and relative abundance r — n/rir (where n is frozen out 
concentration of particles and is concentration of relativistic species) must 
dominate in the Universe before their decay. Dominance of these nonrelativistic 
particles at i > to, where 

to = -TfT' W 
^0 

corresponds to dust like stage with equation of state p = 0, at which particle 
density fluctuations grow as 

S(^t) = ^ cx ^2/3 (7) 

and development of gravitational instability results in formation of gravitation- 
ally bound systems, which decouple at 

t ^ tf ^ USiU)-^^^ (8) 
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from general cosmological expansion, when 6{tf) ^ 1 for fluctuations, entering 
horizon at t = U > to with amplitude d{ti). 

Formation of these systems can result in black hole formation either imme- 
diately after the system decouples from expansion or in result of evolution of 
initially formed nonrelativistic gravitationally bound system. 

If density fluctuation is especially homogeneous and isotropic, it directly 
collapses to BH as soon as the amplitude of fluctuation grows to 1 and the system 
decouples from expansion. A probability for direct BH formation in collapse of 
such homogeneous and isotropic configurations gives minimal estimation of BH 
formation on dust-like stage. 

This probability was calculated in [IT] with the use of the following ar- 
guments. In the period t ~ i/, when fluctuation decouples from expansion, its 
configuration is defined by averaged density pi , size ri , deviation from sphericity 
s and by inhomogeneity u of internal density distribution within the fluctuation. 
Having decoupled from expansion, the configuration contracts and the minimal 
size to which it can contract is 

rmin ~ sri , (9) 

being determined by a deviation from sphericity 

s = max{|7i -72I , I71 - 73| , |73 - 72I}, (10) 

where 71, 72 and 73 define a deformation of configuration along its three main 
orthogonal axes. It was first noticed in [41 that to form a black hole in result 
of such contraction it is sufficient that configuration returns to the size 

Tmin ^ rg ^ U ^ S{U)ri, (11) 

which had the initial fluctuation 6{ti), when it entered horizon at cosmological 
time ti. If 

s < 5{U), (12) 

configuration is sufficiently isotropic to concentrate its mass in the course of 
collapse within its gravitational radius, but such concentration also implies suf- 
ficient homogeneity of configuration. Density gradients can result in gradients 
of pressure, which can prevent collapse to BH. This effect does not take place for 
contracting coUisionless gas of weakly interacting massive particles, but due to 
inhomogeneity of collapse the particles, which have already passed the caustics 
can free stream beyond the gravitational radius, before the whole mass is con- 
centrated within it. Collapse of nearly spherically symmetric dust configuration 
is described by Tolmen solution. It's analysis [321 12041 [55] has provided a 
constraint on the inhomogeneity u = 5pi/ pi within the configuration. It was 
shown that both for coUisionless and interacting particles the condition 

u < S{Uf/^ (13) 

is sufficient for configuration to contract within its gravitational radius. 
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A probability for direct BH formation is then determined by a product of 
probability for sufficient initial sphericity Ws and homogeneity Wu of configu- 
ration, which is determined by the phase space for such configurations. In a 
calculation of Wg one should take into account that the condition ([T^ implies 
5 conditions for independent components of tensor of deformation before its 
diagonalization (2 conditions for three diagonal components to be close to each 
other and 3 conditions for nondiagonal components to be small). Therefore, the 
probability of sufhcient sphericity is given by [HI 12041 [55] 

^ SiU)^ (14) 

and together with the probability for sufhcicnt homogeneity 

Wu 6iuf/^ (15) 

results in the strong power-law suppression of probability for direct BH forma- 
tion 

WpBH = W,-W^^S{uy^/^. (16) 

Though this calculation was originally done in [HI SSI 12041 [55] for Gaussian 
distribution of fluctuations, it does not imply specific form of high amplitude 
tail of this distribution and thus should not change strongly in a case of non- 
Gaussian fluctuations [3T] . 

The mechanism [HI gl IH [20l [35l [3 [4] is effective for formation of PBHs 
with mass in an interval 

Mo < M < Mbh max ■ (17) 

The minimal mass corresponds to the mass within cosmological horizon in the 
period t ^ to, when particles start to dominate in the Universe and it is equal 

to [HHHllllini 133 13 H] 

Mo^^ptl^m.ii'^r. (18) 
3 rm 

The maximal mass is indirectly determined by the condition 

T = t{Mbhmaoc)S{Mbhmaxy^^'^ (19) 

that fluctuation in the considered scale Mbhmax, entering the horizon at t{Mbhmax) 
with an amplitude S{Mbhmax) can manage to grow up to nonlinear stage, de- 
couple and collapse before particles decay at i = r. For scale invariant spectrum 
S{M) = 60 the maximal mass is given by |65| 

Mbhmax = nipi-^Sl^^ = mp;T(5o^^. (20) 
tpi 

The probability, given by Eq. pB]) . is also appropriate for formation of PBHs 
on dust-like preheating stage after inflation [44l [3l |4] . The simplest example 
of such stage can be given with the use of a model of homogeneous massive 
scalar field [3 H] . Slow rolling of the field in the period i <C 1/m (where m is 
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the mass of field) provides chaotic inflation scenario, while at i > l/m the field 
oscillates with period l/m. Coherent oscillations of the field correspond to an 
averaged over period of oscillations dust-like equation of state p = 0, at which 
gravitational instability can develop. The minimal mass in this case corresponds 
to the Jeans mass of scalar field, while the maximal mass is also determined by 
a condition that fluctuation grows and collapses before the scalar field decays 
and reheats the Universe. 

The probability Wpbh{M) determines the fraction of total density 

/3(M) = PEEHiMl ^ WpBH (M) , (21) 

Ptot 

corresponding to PBHs with mass M . For 5{M) ^ 1 this fraction, given by 
Eq. (fTB)l . is small. It means that the bulk of particles do not collapse directly 
in black holes, but form gravitationally bound systems. Evolution of these 
systems can give much larger amount of PBHs, but it strongly depends on 
particle properties. 

Superweakly interacting particles form gravitationally bound systems of col- 
lisionless gas, which remind modern galaxies with coUisionless gas of stars. Such 
system can finally collapse to black hole, but energy dissipation in it and con- 
sequently its evolution is a relatively slow process [2051 [31 Hj. The evolution 
of these systems is dominantly determined by evaporation of particles, which 
gain velocities, exceeding the parabolic velocity of system. In the case of binary 
collisions the evolution timescale can be roughly estimated [2051 [31 3] as 

N . . 

= h^*^^ ^^^^ 

for gravitationally bound system of N particles, where the free fall time tff for 
system with density p is tjf ^ (47rGp)~^/^. This time scale can be shorter due 
to collective effects in coUisionless gas |206) and be at large iV of the order of 

t,,^N^/Hff. (23) 

However, since the free fall time scale for gravitationally bound systems of col- 
lisionless gas is of the order of cosmological time tf for the period, when these 
systems are formed, even in the latter case the particles should be very long 
living r 3> t/ to form black holes in such slow evolutional process. 

The evolutional time scale is much smaller for gravitationally bound systems 
of superheavy particles, interacting with light relativistic particles and radiation. 
Such systems have analogy with stars, in which evolution time scale is defined 
by energy loss by radiation. An example of such particles give superheavy color 
octet fermions of asymptotically free SU(5) model [45j or magnetic monopoles 
of GUT models. Having decoupled from expansion, frozen out particles and 
antiparticles can annihilate in gravitationally bound systems, but detailed nu- 
merical simulation [46 has shown that annihilation can not prevent collapse of 
the most of mass and the timescale of collapse does not exceed the cosmological 
time of the period, when the systems are formed. 
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3.2 Spikes from phase transitions on inflationary stage 

Scale non-invariant spectrum of fluctuations, in which amphtude of small scale 
fluctuations is enhanced, can be another factor, increasing the probability of 
PBH formation. The simplest functional form of such spectrum is represented 
by a blue spectrum with a power law dispersion 

{5\M))oiM-\ (24) 

with amplitude of fluctuations growing at fc > to small M. The realistic 
account for existence of other scalar fields together with inflaton in the period 
of inflation can give rise to spectra with distinguished scales, determined by 
parameters of considered fields and their interaction. 

In chaotic inflation scenario interaction of a Higgs field (j) with inflaton rj can 
give rise to phase transitions on inflationary stage, if this interaction induces 
positive mass term +^ri'^(jp. When in the course of slow rolling the amplitude 
of inflaton decreases below a certain critical value rj^ = m^/v the mass term in 
Higgs potential 

n</>,^) = -^0^ + ^0^ + y^V (25) 

changes sign and phase transition takes place. Such phase transitions on infla- 
tionary stage lead to the appearance of a characteristic spikes in the spectrum of 
initial density perturbations. These spike-like perturbations re-enter the hori- 
zon during the radiation or dust like era and could in principle collapse to form 
primordial black holes. The possibility of such spikes in chaotic inflation sce- 
nario was first pointed out in 207] and realized in [62] as a mechanism of of 
PBH formation for the model of horizontal unification pPHll^l^l^ . 
For vacuum expectation value of a Higgs field 

{^) = j^v (26) 

and A ~ lO^'' the amplitude S of spike in spectrum of density fluctuations, 
generated in phase transition on inflationary stage is given by [62] 

S«± (27) 

with 



where k ~ 1. 

If phase transition takes place at e-folding N before the end of inflation and 
the spike re-enters horizon on radiation dominance (RD) stage, it forms Black 
hole of mass 

M«^exp{2iV}, (29) 
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where Hq is the Hubble constant in the period of inflation. 

If the spike re-enters horizon on matter dominance (MD) stage it should 
form black holes of mass 

2 

Mw^exp{3iV}. (30) 
-no 

4 First order phase transitions as a source of 
black holes in the early Universe 

First order phase transition go through bubble nucleation. Remind the common 
example of boiling water. The simplest way to describe first order phase transi- 
tions with bubble creation in early Universe is based on a scalar field theory with 
two non degenerated vacuum states. Being stable at a classical level, the false 
vacuum state decays due to quantum effects, leading to a nucleation of bubbles 
of true vacuum and their subsequent expansion [212j . The potential energy of 
the false vacuum is converted into a kinetic energy of bubble walls thus making 
them highly relativistic in a short time. The bubble expands till it collides with 
another one. As it was shown in [57l I213j a black hole may be created in a 
collision of several bubbles. The probability for collision of two bubbles is much 
higher. The opinion of the BH absence in such processes was based on strict 
conservation of the original 0(2,1) symmetry. As it was shown in [551 [SDl [CT] 
there are ways to break it. Firstly, radiation of scalar waves indicates the en- 
tropy increasing and hence the permanent breaking of the symmetry during the 
bubble collision. Secondly, the vacuum decay due to thermal fiuctuation does 
not possess this symmetry from the beginning. The investigations [59l [601 E] 
have shown that BH can be created as well with a probability of order unity in 
collisions of only two bubbles. It initiates an enormous production of BH that 
leads to essential cosmological consequences discussed below. 

In subsection 14.11 the evolution of the field configuration in the collisions of 
bubbles is discussed. The BH mass distribution is obtained in subsection l4.2l In 
subsection 14.31 cosmological consequences of BH production in bubble collisions 
at the end of infiation are considered. 

4.1 Evolution of field configuration in collisions of vacuum 
bubbles 

Consider a theory where a probability of false vacuum decay equals F and dif- 
ference of energy density between the false and true vacuum outside equals py. 
Initially bubbles are produced at rest however walls of the bubbles quickly in- 
crease their velocity up to the speed of light v ~ c ^ 1 because a conversion of 
the false vacuum energy into its kinetic ones is energetically favorable. 

Let us discuss dynamics of collision of two true vacuum bubbles that have 
been nucleated in points (ri, ii), (r2, ^2) and which are expanding into false 
vacuum. Following papers [57 ll214j let us assume for simplicity that the horizon 
size is much greater than the distance between the bubbles. Just after collision 
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mutual penetration of the walls up to the distance comparable with its width 
is accompanied by a significant potential energy increase [215] . Then the walls 
reflect and accelerate backwards. The space between them is filled by the field 
in the false vacuum state converting the kinetic energy of the wall back to 
the energy of the false vacuum state and slowdown the velocity of the walls. 
Meanwhile the outer area of the false vacuum is absorbed by the outer wall, 
which expands and accelerates outwards. Evidently, there is an instant when the 
central region of the false vacuum is separated. Let us note this false vacuum bag 
(FVB) does not possess spherical symmetry at the moment of its separation from 
outer walls but wall tension restores the symmetry during the first oscillation of 
FVB. As it was shown in |214| . the further evolution of FVB consists of several 
stages: 

1) FVB grows up to the definite size Dm until the kinetic energy of its wall 
becomes zero; 

2) After this moment the false vacuum bag begins to shrink up to a minimal 
size D*\ 

3) Secondary oscillation of the false vacuum bag occurs. 

The process of periodical expansions and contractions leads to energy losses 
of FVB in the form of quanta of scalar field. It has been shown in the [214U216] 
that only several oscillations take place. On the other hand, important note 
is that the secondary oscillations might occur only if the minimal size of the 
FVB would be larger than its gravitational radius, D* > Vg. Then oscillating 
solutions of "quasilumps" can be realized |217j . The opposite case {D* < Vg ) 
leads to a BH creation with the mass about the mass of the FVB. As it was 
shown in [59l [601 El] the probability of BH formation is almost unity in a wide 
range of parameters of theories with first order phase transitions. 

4.2 Gravitational collapse of FVB and BH creation 

Consider following [SSI [SHI [SD [SSI \M more details the conditions of converting 
FVB into BH. The mass M of FVB can be calculated in a framework of a specific 
theory and can be estimated in a coordinate system K' where the colliding 
bubbles are nucleated simultaneously. The radius of each bubble b' in this 
system equals to half of their initial coordinate distance at first moment of 
collision. Apparently the maximum size Dm of the FVB is of the same order as 
the size of the bubble, since this is the only parameter of necessary dimension 
on such a scale: Dm = 2b'C. The parameter C ~ 1 is obtained by numerical 
calculations in the framework of each theory, but its exact numerical value does 
not affect significantly conclusions. 

One can find the mass of FVB that arises at the collision of two bubbles of 
radius: 

M^^iCb'fp, (31) 

This mass is contained in the shrinking area of false vacuum. Suppose for 
estimations that the minimal size of FVB is of order of wall width A. The BH 
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is created if minimal size of FVB is smaller than its gravitational radius. It 
means that at least at the condition 



A < rg = 2GM (32) 

the FVB can be converted into BH (where G is the gravitational constant) . 
As an example consider a simple model with Lagrangian 

L=\ {d,^f - ^ ('J^' - ^If - + I'o) . (33) 

In the thin wall approximation the width of the bubble wall can be expressed 

as A = 2 (^\f\<^^ . Using ([5^ one can easily derive that at least FVB with 
mass 

M > ^=J— (34) 

should be converted into BH of mass M. The last condition is valid only in case 
when FVB is completely contained within the cosmological horizon, namely 
Mh > l/V^^oG where the mass of the cosmological horizon at the moment 
of phase transition is given by Mh = mp;/$Q. Thus for the potential (|33p at 
the condition A > ($o/™p/)^ a- BH is formed. This condition is valid for any 
realistic set of parameters of theory. 

The mass and velocity distribution of FVBs, supposing its mass is large 
enough to satisfy the inequality ([32|) . has been found in [59l |60l [6T] . This distri- 

bution can be written in the terms of dimensionless mass = (^P) ( J : 

r-TO^ = 647r (f ) ^^e^^^Ji^, v), 
J{^i,v) = J^dTe-^\T^^fi[l+^^l + v)]. 

The numerical integration of psp revealed that the distribution is rather 
narrow. For example the number of BH with mass 30 times greater than the 
average one is suppressed by factor 10^. Average value of the non dimensional 
mass is equal to /i = 0.32. It allows to relate the average mass of BH and volume 
containing the BH at the moment of the phase transition: 

(Mbh) = jm'p^. (Vbh) ^ O.OUp, (Vbh) ■ (36) 

4.3 First order phase transitions in the early Universe 

Inflation models ended by a first order phase transition hold a dignified posi- 
tion in the modern cosmology of early Universe (see for example |218| I219| I220|. 
[2111 muni]). The interest to these models is due to, that such models 
are able to generate the observed large-scale voids as remnants of the primor- 
dial bubbles for which the characteristic wavelengths are several tens of Mpc. 
|223l I224j . A detailed analysis of a first order phase transition in the context of 
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extended inflation can be found in [225j . Hereafter we will be interested only 
in a final stage of inflation when the phase transition is completed. Remind 
that a first order phase transition is considered as completed immediately af- 
ter establishing of true vacuum percolation regime. Such regime is established 
approximately when at least one bubble per unit Hubble volume is nucleated. 
Accurate computation [225] shows that first order phase transition is successful 
if the following condition is valid: 

' lend 

Here T is the bubble nucleation rate. In the framework of first order infiation 
models the filling of all space by true vacuum takes place due to bubble collisions, 
nucleated at the final moment of exponential expansion. The collisions between 
such bubbles occur when they have comoving spatial dimension less or equal 
to the effective Hubble horizon H~^j^ at the transition epoch. If we take i/o = 
lOQhKm/ sec /Mpc in 57 = 1 Universe the comoving size of these bubbles is 
approximately 10^^^h~^ Mpc. In the standard approach it believes that such 
bubbles are rapidly thermalized without leaving a trace in the distribution of 
matter and radiation. However, in the previous subsection it has been shown 
that for any realistic parameters of theory, the collision between only two bubble 
leads to BH creation with the probability closely to 100% . The mass of this 
BH is given by (see (IMI) ) 

Mbh = liMbub (38) 

where 71 ~ 10"^ and M^ub is the mass that could be contained in the bubble vol- 
ume at the epoch of collision in the condition of a full thermalization of bubbles. 
The discovered mechanism leads to a new direct possibility of PBH creation at 
the epoch of reheating in first order inflation models. In standard picture PBHs 
are formed in the early Universe if density perturbations are sufficiently large, 
and the probability of PBHs formation from small post- infiation initial pertur- 
bations is suppressed (see Section [3|). Completely different situation takes place 
at final epoch of first order inflation stage; namely collision between bubbles of 
Hubble size in percolation regime leads to copious PBH formation with masses 

Afo=7a& = V^' (39) 

where Mj!°J^ is the mass of Hubble horizon at the end of inflation. According 
to (|36p the initial mass fraction of this PBHs is given by /Sq ~ 7i/e ~ 6 • 10^'^. 
For example, for typical value of Hend ~ 4 • 10~^mp; the initial mass fraction /3 
is contained in PBHs with mass Mq « Ig. 

In general the Hawking evaporation of mini BHs could give rise to a variety 
possible end states. It is generally assumed, that evaporation proceeds until the 
PBH vanishes completely |226j . but there are various arguments against this 
proposal (see e.g. |227| [75l I228| I229j ). If one supposes that BH evaporation 
leaves a stable relic, then it is naturally to assume that it has a mass of order 
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TTirei = kiTLpi, whcrc fc ~ 1 10^. Wc can investigate the consequences of 
PBH forming at the percolation epoch after first order inflation, supposing that 
the stable relic is a result of its evaporation. As it follows from the above 
consideration the PBHs arc preferentially formed with a typical mass Mq at a 
single time ti. Hence the total density p at this time is 



where /3o denotes the fraction of the total density, corresponding to PBHs in 
the period of their formation ti. The evaporation time scale can be written in 
the following form 

™„ = Ar (41) 

where is the number of effective massless degrees of freedom. 

Let us derive the density of PBH relics. There are two distinct possibilities 
to consider. 

The Universe is still radiation dominated (RD) at tbh- This situation will 
be hold if the following condition is valid Pbh{tbh) < Pj{tbh)- It is possible 
to rewrite this condition in terms of Hubble constant at the end of inflation 

Hend r,5/2 -1/2 i n-6 



rUpl 



> Po g* - 10"'' (42) 



Taking the present radiation density fraction of the Universe to be = 2.5 • 
10~^h~^ {h being the Hubble constant in the units of WOkm - s~^Mpc~^), and 
using the standard values for the present time and time when the density of 
matter and radiation become equal, we find the contemporary densities fraction 
of relics 

Qrel « 1026/i-2fc ( ^) (43) 

V rupi ) 

It is easily to see that relics overdose the Universe i^rei >> 1) for any reasonable 
k and i?end > 10~^mj,z. 

The second case takes place if the Universe becomes PBHs dominated at 
period t\ <t'2,< tbh- This situation is realized under the condition pBH{t2) > 
p-y{t2), which can be rewritten in the form 

< 10-^ (44) 

rripi 

The present day relics density fraction takes the form 

nrel^l0^^h-^k(^Y^^ (45) 

V mpi J 

Thus the Universe is not overdosed by relics only if the following condition is 
valid 

^ < 2 • 10-19/i4/3A:-2/3. (46) 
mpi 
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This condition implies that the masses of PBHs created at the end of inflation 
have to be larger than 

Mo>10"<?-/i'4/3-/fc2/3. (47) 

From the other hand there are a number of well-known cosmological and as- 
trophysical hmits [5301 [HHl HSl 1131 HSHSMIMS] which prohibit the creation 
of PBHs in the mass range ([Tf|) with initial fraction of mass density close to 
Po ~ 10-2. 

So one have to conclude that the effect of the false vacuum bag mechanism 
of PBH formation makes impossible the coexistence of stable remnants of PBH 
evaporation with the first order phase transitions at the end of inflation. 



5 Gravitino production by PBH evaporation and 
constraints on the inhomogeneity of the early 
Universe 

Presently there are no observational evidences, proving existence of PBHs. How- 
ever, even the absence of PBHs provides a very sensitive theoretical tool to study 
physics of early Universe. PBHs represent nonrelativistic form of matter and 
their density decreases with scale factor a as oc oc T^, while the total density 
is oc (X in the period of radiation dominance (RD). Being formed within 
horizon, PBH of mass M, can be formed not earlier than at 

t{M) = tpi = — . (48) 

mpi rripi 

If they are formed on RD stage, the smaller are the masses of PBHs, the larger 
becomes their relative contribution to the total density on the modern MD stage. 
Therefore, even the modest constraint for PBHs of mass M on their density 

npBHiM) = Pf^^ (49) 

Pc 

in units of critical density pc = SH"^ / {SttG) from the condition that their con- 
tribution a{M) into the the total density 

a{M)^P^^^^^n,n,iM) (50) 

Ptot 

for Ptot = Pc does not exceed the density of dark matter 

a{M) = npBH{M) < VLdm = 0.23 (51) 
converts into a severe constraint on this contribution 

(3 ^ P^^"^f\'f^ (52) 

Ptot(tf) 
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in the period t/ of their formation. If formed on RD stage at = t{M), 
given by (|48p . which corresponds to the temperature Tf = rUpi ^/mpi /M, PBHs 
contribute into the total density in the end of RD stage at teq, corresponding to 
Teq « leV, by factor a{teq)/a{tf) = Tf/T^q = mpi/T^qy^mpi/M larger, than in 
the period of their formation. The constraint on /3{M), following from Eq. ((5T|) 
is then given by 




/3(Af)=a(M)^J— <0.23^J— . (53) 

rripi y nipi nipi " ■ 

The possibility of PBH evaporation, revealed by S. Hawking [72], strongly 
influences effects of PBHs. In the strong gravitational field near gravitational 
radius rg of PBH quantum effect of creation of particles with momentum p ~ 
1/rg is possible. Due to this effect PBH turns to be a black body source of 
particles with temperature (in the units fi = c = A: = 1) 

The evaporation timescale BH is tbh ^ 1^% (see Ea. dlD) and discussion 
in previous section) and at M < lO"'^'* g is less, than the age of the Universe. 
Such PBHs can not survive to the present time and the magnitude Ea. ([5T|) for 
them should be re-defined and has the meaning of contribution to the total 
density in the moment of PBH evaporation. For PBHs formed on RD stage and 
evaporated on RD stage at i < tf.q the relationship Eq. (|55)) between (5{M) and 
a{M) is given by [2361135] 

/3(M)=a(Af)^. (55) 

The relationship between (3{M) and a{M) has more complicated form, if PBHs 
are formed on early dust-like stages [33] [35] [75] [3] , or such stages take place 
after PBH formation [75] [3]. Relative contribution of PBHs to total density 
does not grow on dust-like stage and the relationship between (3{M) and a{M) 
depends on details of a considered model. Minimal model independent factor 
a{M) / /3{M) follows from the account for enhancement, taking place only during 
RD stage between the first second of expansion and the end of RD stage at teq, 
since radiation dominance in this period is supported by observations of light 
element abundance and spectrum of CMB [43l [35l [TS] [3] . 

Effects of PBH evaporation make astrophysical data much more sensitive 
to existence of PBHs. Constraining the abundance of primordial black holes 
can lead to invaluable information on cosmological processes, particularly as 
they are probably the only viable probe for the power spectrum on very small 
scales which remain far from the Cosmological Microwave Background (CMB) 
and Large Scale Structures (LSS) sensitivity ranges. To date, only PBHs with 
initial masses between ~ 10^ g and ^ 10^® g have led to stringent limits (see 
e.g. |75l [76] [77] [35]) from consideration of the entropy per baryon, the deu- 
terium destruction, the '*He destruction and the cosmic-rays currently emitted 
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by the Hawking process [75] . The existence of hght PBHs should lead to impor- 
tant observable constraints, either through the direct effects of the evaporated 
particles (for initial masses between 10^"* g and 10^^ g) or through the indi- 
rect effects of their interaction with matter and radiation in the early Universe 
(for PBH masses between 10^ g and 10^^ g). In these constraints, the effects 
taken into account are those related with known particles. However, since the 
evaporation products are created by the gravitational field, any quantum with 
a mass lower than the black hole temperature should be emitted, independently 
of the strength of its interaction. This could provide a copious production of 
superweakly interacting particles that cannot not be in equilibrium with the hot 
plasma of the very early Universe. It makes evaporating PBHs a unique source 
of all the species, which can exist in the Universe. 

Following [3l|4l|78l[66] and [171[|172| (but in a different framework and using 
more stringent constraints), limits on the mass fraction of black holes at the time 
of their formation [(3 = ppsH / ptot) were derived in [237j using the production of 
gravitinos during the evaporation process. Depending on whether gravitinos are 
expected to be stable or metastable, the limits are obtained using the require- 
ment that they do not overdose the Universe and that the formation of light 
nuclei by the interactions of ^He nuclei with nonequilibrium flux of D,T,'^He and 
^He does not contradict the observations. This approach is more constraining 
than the usual study of photo-dissociation induced by photons-photinos pairs 
emitted by decaying gravitinos. It opened a new window for the upper limits on 
f3 below 10^ g. The cosmological consequences of the limits, obtained in [237], 
are briefly reviewed in the framework of three different scenarios: a blue power 
spectrum, a step in the power spectrum and first order phase transitions. 

5.1 Limits on the PBH density 

Several constraints on the density of PBHs have been derived in different mass 
ranges assuming the evaporation of only standard model particles : for 10^ g < 
M < 10^^ g the entropy per baryon at nucleosynthesis was used [138] to ob- 
tain P < (10^ g/M), for 10^ g < M < 10" g the production of nn pairs 
at nucleosynthesis was used [133] to obtain /3 < 3 x 10^1^(10^ g/M)^/^ , 
for 10"'^° g < M < 10^^ g deuterium destruction was used 232 to obtain 
/3 < 3 X 10-22(M/10i° g)i/2^ for 10" g < M < 10" g spallation of ^He 
was used to obtain /3 < 3 x 10~^\M/10^ gf/^, for Af w 5 x 10^^ g the 

gamma-rays and cosmic-rays were used 235' '167' to obtain (3 < 10"^*. Slightly 
more stringent limits were obtained in ^239, . leading to (3 < 10"^" for masses 
between 10^ g and 10" g and in [168], leading to /3 < lO^^s for M = 5 x 10" g. 
Gamma-rays and antiprotons were also recently re-analyzed in [156^ and [153j. 
improving a little the previous estimates. Such constraints, related to phenom- 
ena occurring after the nucleosynthesis, apply only for black holes with initial 
masses above ~ 10^ g. Below this value, the only limits are the very weak en- 
tropy constraint (related with the photon-to-baryon ratio) and the constraint, 
assuming stable remnants of black holes forming at the end of the evaporation 
mechanism as described in the previous Section. 
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To derive a limit in the initial mass range nipi < M < 10^^ g, graviti- 
nos emitted by black holes were considered in |237j . Gravitinos are expected 
to be present in all local supersymmetric models, which are regarded as the 
more natural extensions of the standard model of high energy physics (see, e.g., 
[240j for an introductory review). In the framework of minimal Supergravity 
(mSUGRA), the gravitino mass is, by construction, expected to lie around the 
electroweak scale, i.e. in the 100 GeV range. In this case, the gravitino is 
metastable and decays after nucleosynthesis, leading to important modifications 
of the nucleosynthesis paradigm. Instead of using the usual photon-photino de- 
cay channel, the study of [237] relied on the more sensitive gluon-gluino channel. 
Based on |2411 12421 12431 12441 1245j . the antiprotons produced by the fragmen- 
tation of gluons emitted by decaying gravitinos were considered as a source 
of nonequilibrium light nuclei resulting from collisions of those antiprotons on 
equilibrium nuclei. Then, ^Li, ''Li and ^Be nuclei production by the interactions 
of the nonequilibrium nuclear flux with ^He equilibrium nuclei was taken into 
account and compared with data (this approach is supported by several recent 
analysis |2461 1247] which lead to similar results). The resulting Monte-Carlo 
estimates |244| lead to the following constraint on the concentration of graviti- 
nos: ^3/2 < 1-1 X lO^^'^^^a/V^i where TO3/2 is the gravitino mass in GcV. This 
constraint has been successfully used to derive an upper limit on the reheat- 
ing temperature of the order [244]: Tr < 3.8 x 10^ GeV. The consequences of 
this limit on cosmic-rays emitted by PBHs was considered, e.g., in [169j . In 
the approach of [237] this stringent constraint on the gravitino abundance was 
related to the density of PBHs through the direct gravitino emission. The usual 
Hawking formula 72j was used for the number of particles of type i emitted per 
unit of time t and per unit of energy Q. Introducing the temperature defined 
by Eq. dSJ) T = hc^ /{IGTr'^kGM) « {10^^g)/M GeV, taking the relativistic 
approximation for Fj, and integrating over time and energy, the total number 
of quanta of type i can be estimated as: 

^OT ^ 27 X 10^4 rT,, ^ .X 

where T is in GeV, rupi ~ 10~^ g, a; = Q/T, m is the particle mass and asuGRA 
accounts for the number of degrees of freedom through M^dM = —asuGRAdt 
where M is the black hole mass. Once the PBH temperature is higher than 
the gravitino mass, gravitinos will be emitted with a weight related with their 
number of degrees of freedom. Computing the number of emitted gravitinos as a 
function of the PBH initial mass and matching it with the limit on the gravitino 
density imposed by nonequilibrium nucleosynthesis of light elements leads to an 
upper limit on the PBH number density. If PBHs are formed during a radiation 
dominated stage, this limit can easily be converted into an upper limit on (3 
by evaluating the energy density of the radiation at the formation epoch. The 
resulting limit is shown on Fig. [l]and leads to an important improvement over 
previous limits, nearly independently of the gravitino mass in the interesting 
range. This opens a new window on the very small scales in the early Universe. 
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Figure 1: Constraints of [237j on the fraction of the Universe going into PBHs 
(adapted from [751 [7S1 [771 [3S] ) . The two curves obtained with gravitinos emis- 
sion in mSUGRA correspond to 7713/2 = 100 GeV (lower curve in the high mass 
range) and 7713/2 — 1 TeV (upper curve in the high mass range) 
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It is also possible to consider limits arising in Gauge Mediated Susy Breaking 
(GMSB) models [248]. Those alternative scenarios, incorporating a natural 
suppression of the rate of flavor-changing neutral-current due to the low energy 
scale, predict the gravitino to be the Lightest Supersymmetric Particle (LSP). 
The LSP is stable if R-parity is conserved. In this case, the limit was obtained 
|237j by requiring Vl^/2fl < ^^Af,o> by requiring that the current gravitino 
density does not exceed the matter density. It can easily be derived from the 
previous method, by taking into account the dilution of gravitinos in the period 
of PBH evaporation and conservation of gravitino to specific entropy ratio, that 

/3 < r (57) 



N , "'3/2 



where is the total number of gravitinos emitted by a PBH with initial 

mass M, t^q is the end of RD stage and tf — max(i jo^^^, tg^f;;) when a non- 
trivial equation of state for the period of PBH formation is considered, e.g. a 
dust-like phase which ends at tf-nd IS]- The limit ([57]) does not imply thermal 
equilibrium of relativistic plasma in the period before PBH evaporation and is 
valid even for low reheating temperatures provided that the equation of state 
on the preheating stage is close to relativistic. With the present matter density 
^M,o ~ 0.27 1249J this leads to the limit shown in Fig. [2] for TO3/2 = 10 GeV. 
Following ((57|) this limit scales with gravitino mass as c>c iTi~j^. Models of 
gravitino dark matter with ^3/2,0 — ^cdm,o, corresponding to the case of 
equality in the above formula, were recently considered in [250i 1251] . 



5.2 Cosmological consequences 

Upper limits on the fraction of the Universe in primordial black holes can be 
converted into cosmological constraints on models with significant power on 
small scales }237j . 

The easiest way to illustrate the importance of such limits is to consider a 
blue power spectrum and to derive a related upper value on the spectral index n 
of scalar fluctuations {P{k) oc fc"). It has recently been shown by WMAP |[249j 
that the spectrum is nearly of the Harrison-Zel'dovich type, i.e. scale invariant 
with n K, 1. However this measure was obtained for scales between 10*^ and 
10^*^ times larger that those probed by PBHs and it remains very important to 
probe the power available on small scales. The limit on n given in |237| must 
therefore be understood as a way to constrain P{k) at small scales rather than 
a way to measure its derivative at large scales : it is complementary to CMB 
measurements. Using the usual relations between the mass variance at the PBH 
formation time cr nit form) and the same quantity today an (to) |87j . 

^ (, ) -r (f .( MH{to) \^ ( Mniteq) \^ 
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Figure 2: Constraints of [237j on the fraction of the Universe going into PBHs 
(adapted from |75U76l[77ll35| ). The curve obtained with gravitinos emission in 
GMSB correspond to ?7i3/2 = 10 GeV and scales with gravitino mass as oc m^h- 
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where Mnit) is the Hubble mass at time t and teq is the equilibrium time, it is 
possible to set an upper value on (3 which can be expressed as 



<^H{tfor7n) 



(59) 



where (5mm ~ 0.3 is the minimum density contrast required to form a PBH. The 
limit derived in the previous subsection leads to n < 1.20 in the mSUGRA case 
whereas the usually derived limits range between 1.23 and 1.31 [87l [94l [252] . In 
the GMSB case, it remains at the same level for 7713/2 ~ 10 GeV and is slightly 
relaxed for smaller masses of gravitino. This improvement is due to the much 
more important range of masses probed by the method |237j . 

In the standard cosmological paradigm of inflation, the primordial power 
spectrum is expected to be nearly -but not exactly- scale invariant ^253j . The 
sign of the running can, in principle, be either positive or negative. It has been 
recently shown that models with a positive running a^, defined as 



are very promising in the framework of supergravity inflation (see, e.i/., [254] ). 
The analysis |237j strongly limits a positive running, setting the upper bound 
at a tiny value as < 2 x 10~^. This result is more stringent than the upper limit 
obtained through a combined analysis of Lya forest, SDSS and WMAP data 
|255j . —0.013 < as < 0.007, as it deals with scales very far from those probed 
by usual cosmological observations. The order of magnitude of the running 
naturally expected in most models -either inflationary ones (see, e.g.^ |256j ) or 
alternative ones (see, e.g., [257j )- being of a few times 10"'^ our upper bound 
should help to distinguish between different scenarios. 

In the case of an early dust-like stage in the cosmological evolution [ITJ 
[35l [3l 11] , the PBH formation probability is increased to /3 > J^'^/^ where b is 
the density contrast for the considered small scales (see subsection 13. ip . The 
associated limit on n is strengthened to n < 1.19. 

Following 87J, it is also interesting to consider primordial density perturba- 
tion spectra with both a tilt and a step. Such a feature can arise from underlying 
physical processes |258| and allows investigation of a wider class of inflaton po- 
tentials. If the amplitude of the step is defined so that the power on small scales 
is times higher than the power on large scales, the maximum allowed value 
for the spectral index can be computed as a function of p. Figure [H taken from 
|237j . shows those limits, which become extremely stringent when p is small 
enough, for both the radiation-dominated and the dust-like cases. 

Another important consequence of limits [237j concerns PBH relics dark 
matter (see also discussion in subsection 14. 3|) . The idea, introduced in [259], 
that relics possibly formed at the end of the evaporation process could account 
for the cold dark matter has been extensively studied. The amplitude of the 
power boost required on small scales has been derived, 6.17., in j l73j as a function 




(60) 
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upper limit on n for spectra with a step 
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Figure 3: Upper limit from [237j on the spectral index of the power spectrum 
as a function of the amplitude of the step. 
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of the relic mass and of the expected density. The main point was that the "step" 
(or whatever structure in the power spectrum) should occur at low masses to 
avoid the constraints available between 10^ g and 10^^ g. The limit on (3 derived 
in |237j closes this dark matter issue except within a small window below 10'^ g. 

This result can be re-formulated in a more general way. If the nature of cos- 
mological dark matter is related with superweakly interacting particles, which 
can not be present in equilibrium in early Universe and for which nonequilib- 
rium processes of production e.g. in reheating are suppressed, the early Universe 
should be sufficiently homogeneous on small scales to exclude copious creation 
of these species in miniPBH vaporation. 

Finally, the limits |237| also completely exclude the possibility of a copious 
PBH formation process in bubble wall collisions [59l [60l [61] , considered in the 
previous Section. This has important consequences for the related constraints 
on first order phase transitions in the early Universe and on symmetry breaking 
pattern of particle theory. 

6 Massive Primordial Black Holes from collapse 
of closed walls 

A wide class of particle models possesses a symmetry breaking pattern, which 
can be effectively described by pseudo-Nambu-Goldstone (PNG) field and which 
corresponds to formation of unstable topological defect structure in the early 
Universe (see [65] for review and references). The Nambu-Goldstone nature in 
such an effective description reflects the spontaneous breaking of global U(l) 
symmetry, resulting in continuous degeneracy of vacua. The explicit symmetry 
breaking at smaller energy scale changes this continuous degeneracy by discrete 
vacuum degeneracy. The character of formed structures is different for phase 
transitions, taking place on post-inflationary and inflationary stages. 

6.1 Structures from succession of U(l) phase transitions 

At high temperatures such a symmetry breaking pattern implies the succession 
of second order phase transitions. In the flrst transition, continuous degeneracy 
of vacua leads, at scales exceeding the correlation length, to the formation of 
topological defects in the form of a string network; in the second phase transi- 
tion, continuous transitions in space between degenerated vacua form surfaces: 
domain walls surrounded by strings. This last structure is unstable, but, as was 
shown in the example of the invisible axion |20H 12021 1203] , it is reflected in 
the large scale inhomogeneity of distribution of energy density of coherent PNG 
(axion) field oscillations. This energy density is proportional to the initial value 
of phase, which acquires dynamical meaning of amplitude of axion field, when 
axion mass is switched on in the result of the second phase transition. 

The value of phase changes by 2tt around string. This strong nonhomo- 
geneity of phase, leading to corresponding nonhomogeneity of energy density of 
coherent PNG (axion) field oscillations, is usually considered (see e.g. [260[|261] 
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and references therein) only on scales, corresponduing to mean distance between 
strings. This distance is small, being of the order of the scale of cosmological 
horizon in the period, when PNG field oscillations start. However, since the 
nonhomogeneity of phase follows the pattern of axion string network this ar- 
gument misses large scale correlations in the distribution of oscillations' energy 
density. 

Indeed, numerical analysis of string network (see review in [262]) indicates 
that large string loops are strongly suppressed and the fraction of about 80% 
of string length, corresponding to long loops, remains virtually the same in all 
large scales. This property is the other side of the well known scale invariant 
character of string network. Therefore the correlations of energy density should 
persist on large scales, as it was revealed in |201[ I202i 1203] . 

The large scale correlations in topological defects and their imprints in pri- 
mordial inhomogeneities is the indirect effect of inflation, if phase transitions 
take place after reheating of the Universe. Inflation provides in this case equal 
conditions for phase transition, taking place in causally disconnected regions. 

If phase transitions take place on inflational stage new forms of primordial 
large scale correlations appear. The value of phase after the first phase tran- 
sition is inflated over the region corresponding to the period of inflation, while 
fluctuations of this phase change in the course of inflation its initial value within 
the regions of smaller size. Owing to such fluctuations, for the fixed value of 9qo 
in the period of inflation with e-folding N — 60 corresponding to the part of the 
Universe within the modern cosmological horizon, strong deviations from this 
value appear at smaller scales, corresponding to later periods of inflation with 

< 60. If Oqo < TT, the fluctuations can move the value of Sat to 6* at > tt in 
some regions of the Universe. After reheating in the result of the second phase 
transition these regions correspond to vacuum with O^ac — 27r, being surrounded 
by the bulk of the volume with vacuum 9yac = 0. As a result massive walls are 
formed at the border between the two vacua. Since regions with 9yac = 27r are 
confined, the domain walls are closed. After their size equals the horizon, closed 
walls can collapse into black holes. 

This mechanism can lead to formation of primordial black holes of a whatever 
large mass (up to the mass of AGNs |263| . see for latest review |264j ). Such 
black holes appear in the form of primordial black hole clusters, exhibiting 
fractal distribution in space [267( 12681 165] . It can shed new light on the problem 
of galaxy formation [65l I265( I266j . 

6.2 Formation of closed walls in inflationary Universe 

To describe a mechanism for the appearance of massive walls of a size essentially 
greater than the horizon at the end of inflation, let us consider a complex scalar 
field with the potential [2631 [267l [268l [65] 

V{^) = X{\^\'-f/2r + SV{9), (61) 
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where ip — re . This field coexists with an inflaton field which drives the Hubble 
constant H during the inflational stage. The term 



reflecting the contribution of instanton effects to the Lagrangian renormalization 
(see for example [269 ), is negligible on the inflational stage and during some 
period in the FRW expansion. The omitted term ((62)) becomes significant, when 
temperature falls down the values T ~ A. The mass of radial field component 
r is assumed to be sufficiently large with respect to H, which means that the 
complex field is in the ground state even before the end of inflation. Since the 
term (|62p is negligible during inflation, the field has the form ip « //\/2 • e'^, 
the quantity fO acquiring the meaning of a massless field. 

At the same time, the well established behavior of quantum field fluctua- 
tions on the de Sitter background [18] implies that the wavelength of a vacuum 
fluctuation of every scalar field grows exponentially, having a fixed amplitude. 
Namely, when the wavelength of a particular fiuctuation, in the inflating Uni- 
verse, becomes greater than ^ the average amplitude of this fluctuation 
freezes out at some non-zero value because of the large friction term in the 
equation of motion of the scalar field, whereas its wavelength grows exponen- 
tially. Such a frozen fiuctuation is equivalent to the appearance of a classical 
field that does not vanish after averaging over macroscopic space intervals. Be- 
cause the vacuum must contain fluctuations of every wavelength, inflation leads 
to the creation of more and more new regions containing a classical field of 
different amplitudes with scale greater than H^^ . In the case of an effectively 
massless Nambu-Goldstone field considered here, the averaged amplitude of 
phase fiuctuations generated during each e-fold (time interval is given by 



Let us assume that the part of the Universe observed inside the contemporary 
horizon H^^ = 3000ft.-^Mpc was inflating, over Nu ~ 60 e-folds, out of a single 
causally connected domain of size which contains some average value of 

phase 6*0 over it. When inflation begins in this region, after one e-fold, the 
volume of the Universe increases by a factor e'^ . The typical wavelength of the 
fluctuation 50 generated during every e-fold is equal to . Thus, the whole 
domain containing ^o, after the first e-fold effectively becomes divided into 
separate, causally disconnected domains of size H^^ . Each domain contains 
almost homogeneous phase value 9q ± 59. Thereby, more and more domains 
appear with time, in which the phase differs significantly from the initial value 
^0- A principally important point is the appearance of domains with phase 
9 > n. Appearing only after a certain period of time during which the Universe 
exhibited exponential expansion, these domains turn out to be surrounded by 
a space with phase 6 < tt. The coexistence of domains with phases 9 < tt 
and 6 > TT leads, in the following, to formation of a large-scale structure of 
topological defects. 

The potential ([6T|) possesses a. U{1) symmetry, which is spontaneously bro- 
ken, at least, after some period of infiation. Note that the phase fiuctuations 



5V{9) = A'*(l-cos6'), 



(62) 



59 = H/2TTf. 



(63) 
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during the first e-folds may, generally speaking, transform eventually into fluc- 
tuations of the cosmic microwave radiation, which will lead to imposing restric- 
tions on the scaling parameter /. This difficulty can be avoided by taking into 
account the interaction of the field with the inflaton field (i.e. by making 
parameter / a variable [65j). This spontaneous breakdown is holding by the 
condition on the radial mass, rrir = VXf > H. At the same time the condition 

me = ^ < -ff (64) 

on the angular mass provides the freezing out of the phase distribution until 
some moment of the FRW epoch. After the violation of condition (I64p the 
term (|62|) contributes significantly to the potential (|6ip and explicitly breaks 
the continuous symmetry along the angular direction. Thus, potential (|6ip 
eventually has a number of discrete degenerate minima in the angular direction 
at the points 6„iin = 0, ±27r, ±47r, .... 

As soon as the angular mass mg is of the order of the Hubble rate, the phase 
starts oscillating about the potential minimum, initial values being different 
in various space domains. Moreover, in the domains with the initial phase 
TT < 6 < 2tt, the oscillations proceed around the potential minimum at 6'„ij„ = 
27r, whereas the phase in the surrounding space tends to a minimum at the 
point Omin = 0. Upon ceasing of the decaying phase oscillations, the system 
contains domains characterized by the phase Omin = 27r surrounded by space 
with 9min = 0. Apparently, on moving in any direction from inside to outside of 
the domain, we will unavoidably pass through a point where 6 = tt because the 
phase varies continuously. This implies that a closed surface characterized by 
the phase ^^aH = must exist. The size of this surface depends on the moment 
of domain formation in the inflation period, while the shape of the surface may 
be arbitrary. The principal point for the subsequent considerations is that the 
surface is closed. After reheating of the Universe, the evolution of domains with 
the phase > t: proceeds on the background of the Friedman expansion and is 
described by the relativistic equation of state. When the temperature falls down 
to T* ~ A, an equilibrium state between the "vacuum" phase 9yac = 27r inside 
the domain and the Oyac = phase outside it is established. Since the equation 
of motion corresponding to potential (|62p admits a kink-like solution (see [262] 
and references therein), which interpolates between two adjacent vacua Oyac = 
and Ovac = 27r, a closed wall corresponding to the transition region at = tt is 
formed. The surface energy density of a wall of width ^ 1/rn ^ //A^ is of the 
order of /A^ 0. 

Note that if the coherent phase oscillations do not decay for a long time, 
their energy density can play the role of CDM. This is the case, for example, in 
the cosmology of the invisible axion (see [260|, 126 1' and references therein) . 

It is clear that immediately after the end of inflation, the size of domains 
which contains a phase 9vac > 27r essentially exceeds the horizon size. This 

^ The existence of such domain walls in theory of the invisible axion was first pointed out 
in [270l . 
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situation is replicated in the size distribution of vacuum walls, which appear at 
the temperature T* ~ A whence the angular mass mg starts to build up. Those 
walls, which are larger than the cosmological horizon, still follow the general 
FRW expansion until the moment when they get causally connected as a whole; 
this happens as soon as the size of a wall becomes equal to the horizon size 
Rh- Evidently, internal stresses developed in the wall after crossing the horizon 
initiate processes tending to minimize the wall surface. This implies that the 
wall tends, first, to acquire a spherical shape and, second, to contract toward 
the centre. For simplicity, we will consider below the motion of closed spherical 



The wall energy is proportional to its area at the instant of crossing the hori- 
zon. At the moment of maximum contraction, this energy is almost completely 
converted into kinetic energy [273] . Should the wall at the same moment be 
localized within the gravitational radius, a PBH is formed. 

Detailed consideration of BH formation was performed in [263] . The results 
of these calculations are sensitive to changes in the parameter A and the initial 
phase 9u- As the A value decreases to « IGeV, still greater PBHs appear with 
masses of up to ~ lO""^ g. A change in the initial phase leads to sharp variations 
in the total number of black holes. As was shown above, each domain generates 
a family of subdomains in the close vicinity. The total mass of such a cluster 
is only 1.5-2 times that of the largest initial black hole in this space region. 
Thus, the calculations confirm the possibility of formation of clusters of massive 
PBHs ( ^ IOOM0 and above) in the pregalactic stages of the evolution of the 
Universe. These clusters represent stable energy density fluctuations around 
which increased baryonic (and cold dark matter) density may concentrate in 
the subsequent stages, followed by the evolution into galaxies. 

It should be noted that additional energy density is supplied by closed walls 
of small sizes. Indeed, because the smallness of their gravitational radius, they 
do not collapse into BHs. After several oscillations such walls disappear, leaving 
coherent fluctuations of the PNG field. These fluctuations contribute to a local 
energy density excess, thus facilitating the formation of galaxies. 

The mass range of formed BHs is constrained by fundamental parameters of 
the model / and A. The maximal BH mass is determined by the condition that 
the wall does not dominate locally before it enters the cosmological horizon. 
Otherwise, local wall dominance leads to a superluminal a oa expansion for 
the corresponding region, separating it from the other part of the Universe. This 
condition corresponds to the mass |65j 



The minimal mass follows from the condition that the gravitational radius of 
BH exceeds the width of wall and it is equal to [2671 [55] 



walls 0. 



M, 




(65) 



(66) 



^Tho motion of closed vacuum walls has been driven analytically in [2711 1272] . 
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Closed wall collapse leads to primordial GW spectrum, peaked at 



z/o = 3-10"(A//)Hz 



(67) 



with energy density up to 



Qgw « 10-\f/mpi). 



(68) 



At / lO^^GeV this primordial gravitational wave background can reach 
^Gw ~ 10~®. For the physically reasonable values of 



Another profound signature of the considered scenario are gravitational wave 
signals from merging of BHs in PBH cluster. These effects can provide test of 
the considered approach in LISA experiment. 

7 Discussion 

For long time scenarios with Primordial Black holes belonged dominantly to cos- 
mological anti- Utopias, to "fantasies", which provided restrictions on physics of 
very early Universe from contradiction of their predictions with observational 
data. Even this " negative" type of information makes PBHs an important the- 
oretical tool. Being formed in the very early Universe as initially nonrelativistic 
form of matter, PBHs should have increased their contribution to the total den- 
sity during RD stage of expansion, while effect of PBH evaporation should have 
strongly increased the sensitivity of astrophysical data to their presence. It links 
astrophysical constraints on hypothetical sources of cosmic rays or gamma back- 
ground, on hypothetical factors, causing influence on light element abundance 
and spectrum of CMB, to restrictions on superheavy particles in early Universe 
and on first and second order phase transitions, thus making a sensitive astro- 
physical probe to particle symmetry structure and pattern of its breaking at 
superhigh energy scales. 

Gravitational mechanism of particle creation in PBH evaporation makes 
evaporating PBH an unique source of any species of particles, which can exist 
in our space-time. At least theoretically, PBHs can be treated as source of such 
particles, which are strongly suppressed in any other astrophysical mechanism 
of particle production, cither due to a very large mass of these species, or owing 
to their superweak interaction with ordinary matter. 

By construction astrophysical constraint excludes effect, predicted to be 
larger, than observed. At the edge such constraint converts into an alterna- 
tive mechanism for the observed phenomenon. At some fixed values of param- 
eters, PBH spectrum can play a positive role and shed new light on the old 
astrophysical problems. 



1< A < lO^GeV 



(69) 



the maximum of spectrum corresponds to 

3 • 10"^ < J/Q < 3 • lO^Hz. 



(70) 
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The common sense is to think that PBHs should have small sub-stellar mass. 
Formation of PBHs within cosmological horizon, which was very small in very 
early Universe, seem to argue for this viewpoint. However, phase transitions on 
inflationary stage can provide spikes in spectrum of fluctuations at any scale, or 
provide formation of closed massive domain walls of any size. 

In the latter case primordial clouds of massive black holes around interme- 
diate mass or supermassive black hole is possible. Such clouds have a fractal 
spatial distribution. A development of this approach gives ground for a prin- 
cipally new scenario of the galaxy formation in the model of the Big Bang 
Universe. Traditionally, Big Bang model assumes a homogeneous distribution 
of matter on all scales, whereas the appearance of observed inhomogeneities 
is related to the growth of small initial density perturbations. However, the 
analysis of the cosmological consequences of the particle theory indicates the 
possible existence of strongly inhomogeneous primordial structures in the dis- 
tribution of both the dark matter and baryons. These primordial structures 
represent a new factor in galaxy formation theory. Topological defects such as 
the cosmological walls and filaments, primordial black holes, archioles in the 
models of axionic CDM, and essentially inhomogeneous baryosynthesis (lead- 
ing to the formation of antimatter domains in the baryon-asymmetric Universe 

P71 irm \Tm m7\ irm irm hhoi HSU UHl UHl UHl IMSi [isi n eh H) offer 

by no means a complete list of possible primary inhomogeneities inferred from 
the existing elementary particle models. 

Observational cosmology offers strong evidences favoring the existence of 
processes, determined by new physics, and the experimental physics approaches 
to their investigation. Cosmoparticle physics [Il[ll[3l|4j, studying the physical, 
astrophysical and cosmological impact of new laws of Nature, explores the new 
forms of matter and their physical properties. Its development offers the great 
challenge for theoretical and experimental research. Physics of Primordial Black 
holes can play important role in this process. 
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